Microlensing by natural wormholes: theory and simulations 
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We provide an in depth study of the theoretical peculiarities that arise in effective negative mass 
lensing, both for the case of a point mass lens and source, and for extended source situations. We 
describe novel observational signatures arising in the case of a source lensed by a negative mass. We 
show that a negative mass lens produces total or partial eclipse of the source in the umbra region and 
also show that the usual Shapiro time delay is replaced with an equivalent time gain. We describe 
these features both theoretically, as well as through numerical simulations. We provide negative mass 
microlensing simulations for various intensity profiles and discuss the differences between them. The 
light curves for microlensing events are presented and contrasted with those due to lensing produced 
by normal matter. Presence or absence of these features in the observed microlensing events can 
shed light on the existence of natural wormholes in the Universe. 
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I. INTRODUCTION 

Wormholc solutions to the Einstein field equations 
have been extensively studied in the last decade (see Refs. 
jl],0 and references cited therein, as well as the book by 
Visser ||). Wormholes basically represent bridges be- 
tween otherwise separated regions of the space-time (see 
Fig. 1) and need a special kind of matter in order to ex- 
ist. This matter, known as exotic, violates the energy 
conditions (EC), particularly the null (or averaged null) 
one ||,0 . 

To specify what we are referring to when talking about 
the energy conditions, we shall provide their point- wise 
form. Apart from the null (NEC), they are the weak 
(WEC), the strong (SEC), and the dominant (DEC) 
energy conditions. For a Friedman-Robertson- Walker 
space-time and a diagonal stress-energy tensor = 
(p, — p, — p, —p) with p the energy density and p the pres- 
sure of the fluid, they read: 

NEC (p + p>0), 

WEC ^ 0>0) and (p+p>0), 

SEC ^ (p + 3p> 0) and {p + p > 0), 

DEC ^ (p> 0) and (p±p> 0). (1) 

The EC are, then, linear relationships between the en- 
ergy density and the pressure of the matter generating 
the space-time curvature. We can immediately see why 
the possible violations of the EC are so polemical. If 
NEC is violated, then WEC is also violated. Negative 
energy densities — and so negative masses — are thus phys- 
ically admitted. Nevertheless, it is important to keep in 
mind that the EC of classical General Relativity are only 
conjectures. They are widely used to prove theorems 
concerning singularities and black hole thermodynamics, 



such as the area increase theorem, the topological cen- 
sorship theorem, and the singularity theorem of stellar 
collapse ||. However, all EC lack a rigorous proof and, 
indeed, several situations in which the EC are violated 
are known; perhaps the most quoted being the Casimir 
effect, see Refs. Typically, observed violations are 

produced by small quantum systems, resulting of the or- 
der of h. It is currently far from clear whether there could 
be macroscopic quantities of such an exotic, e.g. WEC- 
violating, matter. If it does exist, macroscopic negative 
masses could be part of the ontology of the universe. 

In fact, the possible existence of negative gravitational 
masses is being investigated at least since the end of the 
nineteenth century . The empirical absence of negative 
masses in the Earth neighborhood could be explained as 
the result of the plausible assumption that, repelled by 
the positive masses prevalent in our region of space, the 
negative ones have been driven away to extragalactic dis- 
tances. Bondi already remarked in |1(J that it is just an 
empirical fact that inertial and gravitational masses are 
both positive quantities. Clearly, no other way better 
than devising observational tests for deciding the con- 
troversy on negative mass existence is available. For in- 
stance, if natural wormholes actually exist in the uni- 
verse (e.g. if the original topology after the Big-Bang 
was multiply connected), then there could be some ob- 
servable electromagnetic signatures that might lead to 
their identification. 

The idea that wormholes can act as gravitational lenses 
and induce a microlensing signature on a background 
source was first suggested by Kim and Sung ]llf . Un- 
fortunately, their geometry was of a perfect alignment of 
a source, both wormhole's mouths and an observer, which 
is, on a common sense ground, quite unlikely. They also 
considered both mouths to be of positive mass. Cramer 
et al. U% carried out more detailed analysis of a negative 
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FIG. 1. Left: Embedding diagram for a wormhole. Two 
mouths, joined by a tunnel, can connect regions otherwise 
very much separated (here the normal space should fold as 
a sheet of paper, whereas the wormhole would be a tunnel 
from one side of the sheet to the other). Right: Embedding 
diagram for a black hole. The singularity here is represented 
as a pinch off of the wormhole tunnel. 

mass wormholes and considered the effect they can pro- 
duce on background point sources, at non-cosmological 
distances. The generalization to a cosmological scenario 
was carried out by Torres et al. although lensing of 
point sources was still used. As far as we are aware, the 
first and only bound on the possible existence of nega- 
tive masses, imposed using astrophysical databases, was 
given by Torres et al. . These authors showed that the 
effective gravitational repulsion of light rays creates two 
bursts, which are individually asymmetric under time re- 
versal, although mirror images of each other. Recently, 
Anchordoqui et al. |l4| searched in existent gamma-ray 
bursts databases for signatures of wormhole microlensing. 
Although they detected some interesting candidates, no 
conclusive results could be obtained. Peculiarly asym- 
metric gamma-ray bursts p"5| , although highly uncom- 
mon, might be probably explained by more conventional 
hypothesis, like precessing fireballs (see, for instance, Ref. 

The case of macrolensing was recently analyzed by us 
as well, showing that if large concentrations of localized 
negative masses do exist, we should be able to see dis- 
tinctive effects upon a deep background jl7|]. All in all, 
possible existence of wormholes or of any kind of nega- 
tive masses in the universe could not yet be discarded. In 
order to confirm or to deny the existence of such masses, 
we need to develop a strong theoretical framework and 
to provide a clear test through observations. 

In this work we revisit the microlensing by natural 
wormholes of stellar and sub-stellar masses. We pro- 
vide an in-depth study of the theoretical peculiarities 
that arise in negative mass microlensing, both for a point 
mass lens and source, and for extended source situations. 
For the first time, we present negative mass microlensing 
simulations, showing the resulting shapes of the images, 
the intensity profiles, the time gain function, the radial 
and tangential magnifications, and other features. In this 
regard, this work extends and deepens previous papers 
in several ways, and gives a complete description from 
where to analyze, from a computational and quantita- 



tive point of view, observational predictions, as the ones 
presented for chromaticity in Ref. |l8| . 

This work is divided as follows. In Section II we in- 
troduce the basic theoretical framework of gravitational 
lensing produced by a generic negative mass, assumed 
here to be a wormhole. Two things should be clear to the 
reader. The first is that the use of wormholes is just to fix 
an interesting theoretical background. Any strut of neg- 
ative mass would produce the same effects. The second 
is that we never consider light going through the worm- 
hole mouths, but being deflected in the neighborhood. 
We introduce the effective refractive index, magnifica- 
tion results, the time gain function, and other features, 
in several subsections of Section II. Section III presents 
microlensing simulations, showing the form and position 
of the produced images. In Section IV, we treat the ex- 
tended source case, also using a numerical code and tak- 
ing into account different source profiles. We finally give 
some concluding remarks. 

II. LENSING BY A POINT NEGATIVE MASS 

In this paper we consider lensing only by a point neg- 
ative mass lens, and thus we can use all the assumptions 
concurrent with the treatment of the Schwarzschild lens: 

• Geometrical optics approximation — the scale over 
which the gravitational field changes is much larger 
than the wavelength of the light being deflected. 

• Small- angle approximation — the total deflection 
angle is small. The typical bending angles involved 
in gravitational lensing of cosmological interest are 
less than < 1'; therefore we describe the lens optics 
in the paraxial approximation. 

• Geometrically-thin lens approximation — the max- 
imum deviation of the ray is small compared to 
the length scale on which the gravitational field 
changes. Although the scattering takes place con- 
tinuously over the trajectory of the photon, the ap- 
preciable bending occurs only within a distance of 
the order of the impact parameter. 

We begin the discussion of gravitational lensing by 
defining two planes, the source and the lens plane. 
These planes, described by Cartesian coordinate systems 
(£1, £2) and (771 , 772), respectively, pass through the source 
and deflecting mass and are perpendicular to the optical 
axis (the straight line extended from the source plane 
through the deflecting mass to the observer). Since the 
components of the image position and the source posi- 
tions are much smaller in comparison to the distances to 
lens and source planes, we can write the coordinates in 
terms of the observed angles. Therefore, the image coor- 
dinates can be written as (61, O2) and those of the source 
as (0i,02)- 



2 



A. Effective refractive index of the gravitational 
field of a negative mass and the deflection angle 

The 'Newtonian' potential of a negative point mass 
lens is given by 



G\M\ 



(& 2 +Z 2)l/2 



(2) 



where b is the impact parameter of the unperturbed light 
ray and z is the distance along the unperturbed light 
ray from the point of closest approach. We have used 
the term Newtonian in quotation marks since it is, in 
principle, different from the usual one. Here the poten- 
tial is positive defined and approaching zero at infinity 
H. In view of the assumptions stated above, we can 
describe light propagation close to the lens in a locally 
Minkowskian space-time perturbed by the positive grav- 
itational potential of the lens to first post-Newtonian or- 
der. In this weak field limit, we describe the metric of the 
negative mass body in orthonormal coordinates a; = ct, 
x = (x l ) by 



ds 2 



I 
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(3) 



where dl = |x| denotes Euclidean arc length. The ef- 
fect of the space-time curvature on the propagation of 
light can be expressed in terms of an effective index of 
refraction n c g [|l9| , which is given by 



n oS = 1 



(4) 



Thus, the effective speed of light in the field of a negative 
mass is 



2 

v cS = c/n cS c + -3> . 

c 



(5) 



Because of the increase in the effective speed of light 
in the gravitational field of a negative mass, light rays 
would arrive faster than those following a similar path 
in vacuum. This leads to a very interesting effect when 
compared with the propagation of a light signal in the 
gravitational field of a positive mass. In that case, light 
rays are delayed relative to propagation in vacuum — the 
well known Shapiro time delay. In the case of a negative 
mass lensing, this effect is replaced by a new one, which 
we shall call time gain. We will describe this effect in 
more detail in the following subsections. 

Defining the deflection angle as the difference of the 
initial and final ray direction 



a 



^OUt ) 



(6) 



where e := dx./dl is the unit tangent vector of a ray x(Z), 
we obtain the deflection angle as the integral along the 
light path of the gradient of the gravitational potential 



a = V±$dl , 



(7) 



where denotes the projection of V<& onto the plane 

orthogonal to the direction e of the ray. We find 



V X *(M) = - 



G\M\h 
(fo 2 + z 2 ) 3 / 2 ' 



Eq. 7 then yields the deflection angle 

AG\M\h 



(8) 



(9) 



It is interesting to point out that in the case of the neg- 
ative mass lensing, the term 'deflection' has its rightful 
meaning — the light is deflected away from the mass, un- 
like in the positive mass lensing, where it is bent towards 
the mass. 



B. Lensing geometry and lens equation 




FIG. 2. Lensing geometry of a negative mass. O is the 
observer, S is the source, W is the negative mass lens, 7i is 
one of the images. (5 is the angle between the source and the 
lens — position of the source, 9 is the angle between the source 
and the image — position of the image, and a is the deflection 
angle, b is the impact parameter and Di, D s and Di s are 
angular diameter distances. Other quantities are auxiliary. 

In Fig. 2 we present the lensing geometry for a point- 
like negative mass. From this figure and the definition of 
the deflection angle (Eq. 6), we can obtain the relation 
between the positions of the source and the image. We 
only need to relate the radial distance of the source and 
the image from the center, since due to circular symme- 
try, the azimuthal angle (p is not affected by lensing. This 
gives 
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(J3 - G)D S = -aD h 



(10) 



(11) 



With the deflection angle (Eq. 9), we can write the lens 
equation as 

, = (+ ffl^ = e+ «M^i. (12) 



C. Einstein radius and the formation of images 

A natural angular scale in this problem is given by the 
quantity 



4G|M| As 
c 2 AA 



(13) 



which is called the Einstein angle. In the case of a 
positive point mass lens, this corresponds to the angle 
at which the Einstein ring is formed, happening when 
source, lens and observer are perfectly aligned. As we 
will see later in this section, this does not happen if the 
mass of the lens is negative. There are other differences 
as well. A typical angular separation of images is of or- 
der 2 #e for a positive mass lens. Sources which are closer 
than about #e to the optical axis are significantly magni- 
fied, whereas sources which are located well outside the 
Einstein ring are magnified very little. All this is different 
with a negative mass lens, but nonetheless, the Einstein 
angle remains a useful scale for the description of the var- 
ious regimes in the present case and, therefore, we shall 
use the same nomenclature for its definition. 

The Einstein angle corresponds to the Einstein radius 
in the linear scale (in the lens plane): 



Re — $eA = 



AGM As A 

~ 2 dT 



(14) 



In terms of Einstein angle the lens equation takes the 
form 



(15) 



which can be solved to obtain two solutions for the image 
position 9: 



'1,2 



p± \fi 2 



40% 



(16) 



Unlike in the lensing due to positive masses, we find that 
there are three distinct regimes here and, thus, can clas- 
sify the lensing phenomenon as follows: 



I: (3 < 2#e There is no real solution for the lens 
equation. It means that there are no images when 
the source is inside twice the Einstein angle. 

II: /3 > 20e There are two solutions, corresponding to 
two images both on the same side of the lens and 
between the source and the lens. One is always in- 
side the Einstein angle, the other is always outside 
it. 

Ill: (3 — 2 9e This is a degenerate case, 9 1.2 = #e! 
two images merge at the Einstein angular radius, 
forming the radial arc (see § HE). 



We also obtain two important scales, one is the Ein- 
stein angle (#e) — the angular radius of the radial critical 
curve, the other is twice the Einstein angle (2 6*e) — the 
angular radius of the caustic. Thus, we have two images, 
one is always inside the 9e, one is always outside; and as 
a source approaches the caustic (2 #e) from the positive 
side, two images coming closer and closer together, and 
nearer and the critical curve, thereby brightening. When 
the source crosses the caustic, the two images merge on 
the critical curve (#e) and disappear. 



D. Time Gain and Time-Offset Function 

Following j2fj , we define a scalar potential, tf)(6), which 
is the appropriately scaled projected Newtonian potential 
of the lens, 



As 2 



$(A#, z)dz 



A A c 2 

For a negative point mass lens it is 
AA c l 



(17) 



(18) 



The derivative of ip with respect to 8 is the deflection 
angle 

V e V - DiV b i> = l^i / V^dz = a . (19) 

Thus, the deflection angle is the gradient of ip — the de- 
flection potential 

a(0) = Vg^j . (20) 
^From this fact and the lens equation (11) we obtain 

(0-(3) + V g iP(0) = Q . (21) 
This equation can be written as a gradient, 
1 



= 



(22) 



If we compare this equation with that for the Fermat's 
principle |£0[ 
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V fl i(0) = 0, 



(23) 



wc see that wc can define the time-offset function (oppo- 
site to time-delay function in the case of positive mass 
lens) as 



t{6) 



(1 + zi) AA 
c As 



i(. 



P) 2 + H0) 



+ t 



pot ) 



(24) 



where < gcom is the geometrical time delay due to the ex- 
tra path length of the deflected light ray relative to the 
unperturbed one. It remains the same as in the positive 
case — increase of light-travel-time relative to an unbent 
ray. The coefficient in front of the square brackets en- 
sures that the quantity corresponds to the time offset as 
measured by the observer. The second term t po t is the 
time gain a ray experiences as it traverses the deflection 
potential ip(6), with an extra factor (1 + z{) for the cos- 
mological 'redshifting'. Thus, cosmological geometrical 
time delay is 



(i 



c As 2 



and cosmological potential time gain is 
(1 + zi) AA 



t 



pot 



A, 



(25) 



(26) 



In Fig. |^ we show the time delay and time gain func- 
tions. The top panel shows t ge0 m for a slightly offset 
source. The curve is a parabola centered on the position 
of the source, f3. The central panel displays t po t for a 
point negative mass lens. This curve is centered on the 
lens. The bottom panel shows the total time-offset. Im- 
ages are located at the stationary points of t to tai- Here 
we see two extrema — maximum and minimum — on the 
same side (right) from the optical axis (marked by dots). 

We can find the time difference between the two im- 
ages, 6 1 and #2, that is, if a source has intrinsic variability, 
it will appear in the two images at an interval 



-(1 + zO (V 



/2 



(27) 



where by v we denoted the ratio of absolute values of 
magnifications of images, 



Mi 

M2 



(28) 



and u is the scaled angular position of the source u 
P/6e and r s is the Schwarzschild radius of the lens. 



E. Magnifications 

Light deflection not only changes the direction but also 
the cross-section of a bundle of rays. For an infinites- 
imally small source, the ratio between the solid angles 
gives the flux amplification due to lensing 




pot 




angular position 

FIG. 3. Geometric time delay, gravitational time gain and 
total time offset produced by a point negative mass lens for a 
source that is slightly off the optical axis. 
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du>i 



(29) 



For an infinitesimal source at angular position (3 and im- 
age at angular position 0, the relation between the two 
solid angles is determined by the area distortion, given 
in turn by the determinant of the Jacobian matrix A of 
the lens mapping (3 



A 



d/3 



For a point mass lens magnification is given by 



(3df3 
9d6 



(30) 



(31) 



The image is thus magnified or demagnified by a factor 
\fi\. If a source is mapped into several images, the total 
amplification is given by the sum of the individual image 
magnifications. From the lens equation (15), we find 



Thus, 



e 2 



Ml,2 = 



dp 

dd 



2 



1 - 



n, 2 



(32) 



(33) 



and using u from (28), we find the total magnification 
(Fig. 4, bottom panel, continuous curve) as 
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Mtot = I Mi I + |A*2 1 



u\f u 2 — 4 



(34) 



The tangential and radial critical curves follow from 
the singularities in 



A*tan — 



and 



/^rad 



(35) 



(36) 



/i ra d diverges when = #e — angular radius of the radial 
critical curve. /Li tan always remains finite, which means 
that there are no tangential critical curves — no tangen- 
tial arcs can be formed by the negative point mass lens. 
In Fig. 4 we show the magnification curves (radial, tan- 
gential and total) for both positive (upper panel) and 
negative mass lenses (bottom panel). The difference can 
be seen as follows — in the upper panel there is no singu- 
larity in the radial curve (no radial arcs are formed by the 
positive mass lens), whereas in the bottom panel we see 
that the curve for the radial magnification experiences a 
singularity. 



III. MICROLENSING 



When the angular separation between the images AO 



AO 



A6l 



(37) 



is of the order of milliarcsecs, we cannot resolve the two 
images with existing telescopes and we can only observe 
the lensing effect through their combined light intensity. 
This effect is called microlensing. Both the lens and the 
source are moving with respect to each other (as well as 
the observer). Thus, images change their position and 
brightness. Of particular interest are sudden changes in 
luminosity, which occur when a compact source crosses a 
critical curve. For the positive mass lensing the situation 
is quite simple (Fig. ^ (for a review on the positive mass 
microlensing and its applications, see p6| ). 

For a negative mass lens the situation is different. We 
define a dimensionless minimum impact parameter Bq, 
expressed in terms of the Einstein radius, as the shortest 
distance between the path line of the source and the lens. 
For three different values of Bo we have three different 
lensing configurations shown in Figs. 6, 7, and 8. Note 
the large difference in the shapes of the images for these 
three regimes. In Fig. 9 we show the case of a minimum 
impact parameter equal to zero, B$ = 0, that is, the path 
of the source passes through the lens. 

It can be assumed without the loss of generality that 
the observer and the lens are motionless and the source 
moves in the plane perpendicular to the line of sight 




6 [arcsec] 



9 foresee] 

FIG. 4. The magnifications: tangential /i ta n (dotted lines), 
radial /i ra d (dash-dotted lines), and total fi (continuous 
curves), are plotted as functions of the image position 9 for 
two cases; in the upper panel for the positive mass, in the 
bottom panel for the negative mass. The singularities of fj, ta . n 
and /i r ad give the positions of the tangential and radial critical 
curves, respectively. In the upper panel the singularity is in 
the tangential critical curve, in the bottom panel, instead, in 
the radial critical curve. Here \M\ = 1M , D s = 0.05 Mpc 
and D\ = 0.01 Mpc. Angles are in arcseconds. 
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FIG. 5. Schematic representation of the geometry of the 
positive mass lensing due to the motion of the source, lens and 
the observer (in this case we can consider only the motion of 
the source in the plane perpendicular to the optical axis). The 
lens is indicated with a dot at the center of the Einstein ring, 
which is marked with a dashed line. The positions of the 
source center are shown with a series of small open circles. 
The locations and the shapes of the two images are shown 
with a series of dark ellipses. At any instant, the two images, 
the source and the lens are all on a single line, as shown in 
the figure for one particular instant. 




FIG. 7. True motion of the source and apparent motion of 
the images for Bg = 2. The inner dashed circle is the Einstein 
ring, the outer dashed circle is twice the Einstein ring. The 
rest is as in Fig. S 




FIG. 8. True motion of the source and apparent motion of 
the images for Bo < 2. The inner dashed circle is the Einstein 
FIG. 6. True motion of the source and apparent motion of ring) the outer dashed circle is twice the Einstein ring. The 
the images for Bq > 2. The inner dashed circle is the Einstein res ^ j s as j n pjg |^| 
ring, the outer dashed circle is twice the Einstein ring. The 
rest is as in Fig. pi 
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FIG. 9. True motion of the source and apparent motion 
of the images for Bo — 0. The inner dashed circle is the 
Einstein ring, the outer dashed circle is twice the Einstein 
ring. Images here are shown with the shaded ellipses. The 
rest is as in Fig. ^ 

(therefore, changing its position in the source plane). We 
adopt the treatment given in [^jj for the velocity V, and 
consider effective transverse velocity of the source rela- 
tive to the critical curve (see Appendix A). We define the 
time scale of the microlensing event as the time it takes 
the source to move across the Einstein radius, projected 
onto the source plane, £o = 9eD s , 



t = — 



V 



Angle (3 changes with time according to 




(38) 



(39) 



Here the moment t = corresponds to the smallest an- 
gular distance (3 between the lens and the source. Nor- 
malizing (39) to 0e, 



«(*) = 



' Vt 
0v.D f 



(40) 



where dimensionless impact parameter u is defined in 
(28). Including the time scale t v and defining 



Bo 



0e 



we obtain 



«(*) = A/-Bg+ ( r 



(41) 



(42) 



Finally, the total amplification as a function of time is 
given by 



A(t) = 



(43) 



Comparing this analysis with that of Cramer et al. |L2) , 
we must note that they wrote the equation for the time 
dependent dimensionless impact parameter as (cf. our 
Eq. P) 



B(t) = B \l 



and defined the time scale for the microlensing event as 
the time it takes to cross the minimum impact parameter 
(cf. our Eq. |3§| ) 

^ _ bo 

where bo is the minimum impact parameter and other 
variables carry the same meanings as in our paper. While 
there is no mistake in using such definitions, there is a 
definite disadvantage in doing so. Using Eq. 10 of for 
B(t) we cannot build the light curve for the case of the 
minimum impact parameter Bo = 0. In this case their 
Eq. 8 diverges, although there is nothing wrong with this 
value of Bo (see our Figs. 9 and 10). In the same way, 
their definition of a time scale does not give much infor- 
mation on the light curves. With our definition (Eq. |36| ) 
we can see in Fig. 10 that in the extreme case of Bo = 
the separation between the half-events is exactly 2#e! it 
is always less than that with any other value of Bo- 

In Fig. [To] we show the light curves for the point source 
for four source trajectories with different minimum im- 
pact parameters Bo- As can be seen from the light curves, 
when the distance from the point mass to the source tra- 
jectory is larger than 26b, the light curve is identical to 
that of a positive mass lens light curve. However, when 
the distance is less than 26*e (or in other terms, Bq < 2.0), 
the light curve shows significant differences. Such events 
are characterized by the asymmetrical light curves, which 
occur when a compact source crosses a critical curve. A 
very interesting, eclipse-like, phenomenon occurs here; a 
zero intensity region (disappearance of images) with an 
angular radius 0o 



% = \/40? 



y E A) ' 

or in terms of normalized unit #e , 



A = \ / 4 — £?q 



(44) 



(45) 



In the next Section we shall see how these features get 
affected by the presence of an extended source. 



IV. EXTENDED SOURCE 

In the previous sections we considered magnifications 
and light curves for point sources. However, sources are 
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FIG. 10. Light curves for the negative mass lensing of a 
point source. ^From the center of the graph towards the cor- 
ners the curves correspond to Bo = 2.5, 2.0, 1.5, 0.0. The time 
scale here is £o divided by the effective transverse velocity of 
the source. 



extended, and although their size may be small compared 
to the relevant length scales of a lensing event, this exten- 
sion definitely has an impact on the light curves, as will 
be demonstrated below. From variability arguments, the 
optical and X-ray continuum emitting regions of quasars 
are assumed to be much less than 1 pc |22| , whereas the 
broad- line emission probably has a radius as small as 0.1 
pc [ p3[ . The high energy gamma-spheres have a typical 
radius of 10 15 cm f24| . Hence, one has to consider a fairly 
broad range of source sizes. 

We define here the dimensionless source radius, R, as 



R 



P_ 

0e 



R 

To 



(46) 



where p and R is the angular and the linear physical size 
of the source, respectively, and £o is the the length unit 
in the source plane (see Eq. 38). 



A. Comments on numerical method and simulations 

It is convenient to write the lens equation in the scaled 
scalar form 



where we normalized the coordinates to the Einstein an- 
gle: 



y = 



(48) 



The lens equation can be solved analytically for any 
source position. The amplification factor, and thus the 
total amplification, can be readily calculated for point 
sources. However, as we are interested in extended 
sources, this amplification has to be integrated over the 
source (Eq. fl9]), and furthermore, as we want to build 
the light curves, the total amplification for an extended 
source has to be calculated for many source positions. 
The amplification A of an extended source with surface 
brightness profile I{y) is given by 



A 



J riV(y)A(y) 
JcPyI(y) 



(49) 



where -4o(y) is the amplification of a point source at 
position y. We have used the numerical method first de- 
scribed in |2j|. We cover the lens plane with a uniform 
grid. Each pixel on this grid is mapped, using Eq. H7L 
into the source plane. The step width (5000 x 5000jis 
chosen according the desired accuracy (i.e. the observ- 
able brightness). For a given source position (j/oi,J/02) 
we calculate the squared deviation function (SDF) 



S 2 = (yw - yi(xi,x 2 )) 2 + (2/20 - y2{xi,x 2 )) 2 



(50) 



y = x + 



(47) 



The solutions of the lens equation (Eq. |4g) are given by 
the zeroes of the SDF. Besides, Eq. [30] describes circles 
with radii S around (t/10,2/20) m the source plane. Thus, 
the lines S = const are just the image shapes of a source 
with radius S, which we plot using standard plotting soft- 
ware. Therefore, image points where SDF has value S 2 
correspond to those points of the circular source which 
are at a distance S from the center. The surface bright- 
ness is preserved along the ray and if I(Rq) = lo for the 
source, then the same intensity is given to those pixels 
where SDF= Rq. In this way an intensity profile is cre- 
ated in the image plane and integrating over it we can 
obtain the total intensity of an image. Thus, we obtain 
the approximate value of the total magnification by esti- 
mating the total intensity of all the images and dividing 
it by that of the unlensed source, according to the corre- 
sponding brightness profile of the source (see Appendix 
B). For a source with the constant surface brightness the 
luminosity of the images is proportional to the area en- 
closed by the line S = const. And the total magnification 
is obtained by estimating the total area of all images and 



1 Note, that for the case in which x and y are expressed 
in length units, we obtain a different normalization in each 
plane, which is not always convenient. 
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dividing it by that of the unlensed source. For calcula- 
tions of light curves we used a circular source which is 
displaced along a straight line in the source plane with 
steps equal to 0.01 of the Einstein angular radius. 

B. Results 

In Figs. 11 and 12 we show four projected source 
and image positions, critical curves/caustics in the 
lens/source plane and representative light curves for dif- 
ferent normalized source sizes. The sources are taken to 
be circular disks with constant surface brightness. In or- 
der to get absolute source radii and real light curves we 
need the value of 9e, the normalized angular unit, the 
distance to the source, as well as the velocity V of the 
source relative to the critical curves in the source plane 
(see Appendix A). We have used M = M Q , H Q = 100 
km s _1 Mpc -1 and a standard cosmological model with 
zero cosmological constant. Here and in all subsequent 
simulations the redshift of the source is z s = 0.5 and the 
redshift of the lens is z\ = 0.1. 

We display two cases for two different impact param- 
eters. It must be noted here that the minimum impact 
parameter B defines now the shortest distance between 
the line of path of the center of the source and the lens. 
For each one of B , the dimensionless radius of the source 
R increases from 0.01 to 2.0 in normalized units of #e- 
The shape of produced images changes notably with the 
increase of the source size, as can be seen in the bottom 
right panel of Figs. 11 and 12. At the same time the 
smaller the source the greater the magnification, since 
when the source radius is greater than the Einstein ra- 
dius of the lens, the exterior parts, which arc amplified, 
compete with the interior ones, which are demagnificd. 

It can be noted that, despite the noise in some of the 
simulated light curves, the sharp peaks which occur when 
the source is crossing the critical line are well defined 
even for the smallest source. Note that all infinities are 
replaced now by finite amplifications, and that the curves 
are softened; all these effects being generated by the fi- 
nite size of the source. Indeed, while the impressive drop 
to zero in the light curve is maintained, the divergence to 
infinity, that happens for a point source, is very much re- 
duced. Note, that in cases of a large size of the source the 
magnification is very small. If we would like to see big- 
ger enhancement than that, we should consider sources 
of smaller sizes, approaching the point source situation 
(cf. Fig. 11, upper left plot). 

It is also interesting to note here that for the im- 
pact parameter Bq = 2.0 the light curve of a small ex- 
tended source, though approaching the point source pat- 
tern (Fig. 10), still differs considerably from it (Fig. 12, 
upper left plot). 

In Figs. 13 and 14 we show the images of an extended 
source with a Gaussian brightness distribution for two 
effective dimensionless source radii (Appendix B), Rs, 



-5 5 -5 5 




FIG. 11. Four sets of lens-source configurations (upper pan- 
els) and corresponding amplification as a function of source's 
center position (bottom panels) are shown for four different 
values of the dimensionless source radius R (0.01, 0.5, 1.0, 
2.0, in normalized units, #e). Each of the four upper panels 
display the time dependent position of the source's center, the 
shapes of images (shaded ellipses) and critical curves (dashed 
circles). The series of open small circles show the path of 
the source center. The lens is marked by the central cross. 
Minimum impact parameter Bo — 0. By replacing £ with 
^e-DsV^ 1 = £t v we get corresponding time depending light 
curve. 
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FIG. 12. Same as in Fig. 11, but for minimum impact pa- 
rameter Bq — 2.0. 



of 3.0 and 1.0 (Fig. 13, frames a-e and Fig. 14, frames 
a-e, respectively), together with the corresponding light 
curves (Fig. 13, frame f and Fig. 14, frame f, respec- 
tively). Here the source path passes through the lens 
(Bq = 0), which lies exactly in the center of each frame. 
In Fig. 13 the source's extent in the lens plane is greater 
than the Einstein radius of the lens. Annotated wedges 
provide colour scale for the images. We notice there that 
there is an eclipse-like phenomenon, occurring most no- 
tably when most of the source is near or exactly behind 
the lens. This is consistent with the light curve (frame 
f), where there is a de-magnification. For the source with 
radius smaller than the double Einstein radius of the lens 
(Fig. 14), the low intensity region is replaced by the zero 
intensity region; the source completely disappears from 
the view (frame c). 

In Fig. 15 we display the images of the source with 
the exponential brightness profile (see Appendix B) and 
the corresponding light curve (frame f). The effective 
radius of the source is 1.5. The impact parameter here 
is Bo = 2.0; the lensing regime corresponds to the one 
schematically depicted in Fig. 7. We see how shapes of 
the images change, becoming elongated and forming the 
radial arc (frames c and d). 

In order to compare a constant surface brightness 
source with more realistic distributions, we simulate im- 
ages configurations and calculate light curves for two dif- 
ferent assumed profiles with radial symmetry (see Ap- 
pendix B). 

In Fig. |l6| we compare light curves for three differ- 
ent radially symmetric source profiles, uniform, Gaussian 
and exponential, for two dimensionless source radii R = 
flf" = R c s xpon = 0.1 and R = Rf uss = i?° xpon = 1.0. 
As a reference curve we show the light curve of the point 
source. All curves are made for the impact parameter 
Bq = 0. We can see the larger noise in the uniform 
source curve, since the source with uniform brightness 
has extremely sharp edge, whereas Gaussian and expo- 
nential sources are extremely smooth. Though we con- 
sidered the sources with the same effective radius, we 
can see from the plot that for a small source size, the 
maximum magnification is reached by the source with 
exponential profile (upper panel), which is explained by 
the fact that this profile has a more narrow central peak 
than the Gaussian. 

For the larger source, this behaviour smoothens, 
though we still can see large differences in the light curves 
(bottom panel), where the uniform source experiences 
darkening, while sources with other profiles only undergo 
demagnification. 



V. CONCLUDING REMARKS 



In this paper we have explored the consequences, re- 
garding gravitational microlensing, of the existence of 
matter violating the energy conditions. We have also 
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FIG. 13. Image configurations (frames a to e) and a corre- 
sponding light curve (f) for a Gaussian source with effective 
radius Rs = 3.0, in units of the Einstein angle. The source is 
moving from the lower left corner (frame a) to the right upper 
corner (frame e), passing through the lens (Bo = 0). The lens 
is in the center of each frame. Size of each frame is 5 x 5, 
in the normalized units. Wedges to each frame provide the 
colour scale for the images. Note the eclipse-like phenomenon, 
consistent with the incomplete demagnification showed in the 
light curve (frame f). 



FIG. 14. Image configurations (frames a to e) and a corre- 
sponding light curve (f) for a Gaussian source with effective 
radius Rs = 1.0, in units of the Einstein angle. The source is 
moving from the lower left corner (frame a) to the right upper 
corner (frame e) , passing through the lens (Bo = 0) . The lens 
is in the center of each frame. Size of each frame is 3 x 3, 
in the normalized units. Wedges to each frame provide the 
colour scale for the images. Note the complete disappearence 
of the source when it is inside the double Einstein radius of 
the lens (frame c), corresponding to the drop of magnification 
to zero in the light curve (frame f). 
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FIG. 15. Image configurations (frames a to e) and a cor- 
responding light curve (f) for a source with the exponential 
brightness distribution and effective radius Rs = 1.0, in units 
of the Einstein angle. The source is moving from the upper 
left corner (frame a) to the upper right corner (frame e) with 
the impact parameter Bo = 2.0. The lens is in the center of 
each frame. Size of each frame is 2.5 x 2.5, in the normalized 
units. Wedges to each frame provide the colour scale for the 
images. 





f 



FIG. 16. Light curves for the point mass source (dashed 
line), source with constant surface brightness (solid), source 
with Gaussian brightness distribution (dash-dotted) and ex- 
ponential brightness distribution (dotted) for two different ef- 
fective dimensionless source radii, 0.1 (upper panel) and 1.0 
(bottom panel). 



quantitatively analyzed, using numerical simulations, the 
influence of a finite size of the source on the gravitational 
lensing negative-mass event. We have thus enhanced and 
completed previous works, where the focus was put on 
the point source light curves and no discussion was given 
concerning the shapes of images, actual simulations of 
microlensing events, time gain function, and other fea- 
tures presented here. Figs. 4, 6-9 and 11-16 comprise our 
new results: a useful comparison arena where to test ob- 
servationally the possible existence of wormholes or any 
other kind of negative mass compact objects. 

The next step would be to test these predictions us- 
ing archival, current, and forthcoming observational mi- 
crolensing experiments. The only search done up to now 
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included the BATSE database of gamma-ray bursts, but 
there is still much unexplored territory in the gravita- 
tional microlensing archives. We suggest to adapt the 
alert systems of these experiments in order to include the 
possible effects of negative masses as well. This, perhaps, 
would not imply a burdensome work, but there could be 
a whole new world of discoveries. 
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APPENDIX A: VELOCITIES AND TIME SCALES 

Let the source have a transverse velocity v s measured 
in the source plane, the lens a transverse velocity Vi mea- 
sured in the lens plane, and the observer a transverse ve- 
locity v bs measured in the observer plane. The effective 
transverse velocity of the source relative to the critical 
curves with time measured by the observer is 



1 



1 



1 D s 
l + z\D } 



■vi 



1 As 
l + z\ D] 



■Vobs 



(Al) 



where this effective velocity is such that for a stationary 
observer and lens, the position of the source will change 
in time according to 5£ = VAt. 

We basically have two time scales of interest here. The 
first one is the typical rise time to a peak in the ampli- 
fication. We can estimate that it corresponds to a dis- 
placement of Ay ~ R of the source across a critical line; 
the corresponding time scale n is 



t v R , (A2) 
or in terms of the physical 



where t v is given by (|3 
source size R = £,qR, 



R 

Tl = v 



(A3) 



where effective transverse velocity of the source V is given 
by (Al). The second time scale of interest is the time 
between two peaks r 2 . For a point source we can estimate 
it as 



TABLE I. Time scales for several source radii. R is 
the dimensionless source radius, in units of Einstein angle, 
\M\ = 1.0 M©, redshift of the lens is z\ = 0.1, redshift of the 
source is z s = 0.5, £o = 5.42 x 10 11 km is the normalized 
length unit in the source plane. The time scales correspond 
to apparent source velocity (see Eq. Al) V = 5000 km s _1 . 



R 


R (pc) 


Afi2 a (sec) 


n (yr) 


T2 h (yr) 


0.0 


point source 


2.0 x 10" 4 




6.78 


0.01 


1.07 x 10" 4 




0.03 


6.81 


0.1 


1.07 x 10 -3 




0.34 


7.11 


1.0 


1.07 x 10" 2 




3.38 


10.16 


2.0 


0.3 x 10" 1 




6.75 


13.6 



a u=4.0 (definition in Eq. ||) 
b B = 0.0 



where A is given in Eq. (|4^). For a source with radius R 
and impact parameter Bq = 0.0 it can be shown to be 



T 2 = U 



R 



(A5) 



In Table I we list for different values of source radii the 
time scales t\ and t%, and the time delay between two 
images Atn for the point source. The value of V we 
estimate to be V = 5000 km s _1 . 



APPENDIX B: EXTENDED SOURCES 
BRIGHTNESS PROFILES 

1. Source of uniform brightness 

For a circular source of radius R and uniform bright- 
ness, equation (49) transforms into 



A = 



J d 2 yI(y)A (y) 
nR 2 



(Bl) 



2. Source with Gaussian brightness distribution 



For a Gaussian source we have 
I(r) = I e- r2 / r o , 



(B2) 



where we normalized the profile such that the maximum 
value of / equals unity. We define the radius containing 
90% of all the luminosity as the effective radius of the 
source, Rs- To find the relation between Rg and ro, we 
write the total luminosity as 



L(oo) 



(B3) 



r 2 = t v A 



(A4) 



L(R) 



— r 2 /r 2 r> t 2 

e ' zirr ar = 7rr 



, (B4) 
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then 



L(Rs) 
L(oo) 



= 0.9 = 



1 _ e -«s/'o 



Thus, effective radius relates to the parameter r$ as 
Rs 



= Vln 10 



3. Source with exponential brightness distribution 

We have 

J(r) = I e~ r/ro . (B7) 

In the same manner as above, Rg is defined as radius, 
containing 90% of total luminosity. In the same way as 
above, total luminosity 



L(oo) 



e -r/ro 27rr dr = 2nrl , 



then 
L(R) 

and 



e- r/r °2nrdr = 2tt 



(Rr + r%)e 



L(Rs) 
L(oo) 



= 0.9 . 



From where we find that effective radius relates to the 
parameter r as 



To 



0.1 



Solution to this equation gives i?s/ r o ~ 3.89. This pro- 
file is also normalized such that the maximum value of / 
equals unity. 
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